Spin-ordering and magnetoelastic coupling in the extended Kagome system YBaCo407 
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Low temperature magnetic and structural behavior of the extended Kagome system YBaCo407 
has been studied by single crystal neutron diffraction and high-resolution powder X-ray diffraction. 
Long-range magnetic ordering associated with a structural transition from orthorhombic Pbn2i to 
monoclinic P2i symmetry has been found at Ti ^ 100 K. The interplay between the structural and 
magnetic degrees of freedom testifies that the degeneracy of the magnetic ground state, present in 
the orthorhombic phase, is lifted through a strong magnetoelastic coupling, as observed in other 
frustrated systems. At T2 '^ 60 K, an additional magnetic transition is observed, though iso- 
symmetric. Models for the magnetic structures below Ti and T2 are presented, based on refinements 
using a large number of independent reflections. The results obtained are compared with previous 
single crystal and powder diffraction studies on this and related compositions. 

PACS numbers: 75.25.-j, 75.85.+t, 61.05.F- 
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I. INTRODUCTION 

Magnetic behavior of systems with competitive ex- 
change interactions has been the focus of attention over 
the last several decades. This competition can arise 
due to comparable strength between nearest and beyond 
nearest neighbor interactions or due to specific symme- 
try of the exchange topology.^ The latter, usually re- 
ferred to as geometrical frustration, can be found in 
many materials where magnetic ions form a network of 
triangular-based structural units, as for example in two- 
dimensional Kagome and three-dimensional pyrochlore 
lattices.^ The hallmark of frustration is a large magnetic 
ground state degeneracy and consequently the suppres- 
sion of long-range magnetic ordering. This degeneracy is 
possibly lifted by structural distortions^^ either associ- 
ated to non magnetic instability (bonding requirements, 
orbital/charge order) or as a consequence of strong mag- 
netoelastic coupling. In the resulting distorted phases, 
complex non-collinear or partially disordered spin config- 
urations can be observed. ^^^ Recently, this kind of behav- 
ior has been reported for i?BaCo407 (i?-rare earth or Y) 
cobaltites.^'^ These compositions belong to a new class 
of frustrated systems with trigonal symmetry (P31c), 
formed by an alternate stacking of Kagome and triangu- 
lar layers along the c-axis.^^ The unique exchange topol- 
ogy is formed by a network of corner-sharing trigonal 
bipyramids and triangular clusters ii A structural phase 
transition from trigonal (P31c) to orthorhombic symme- 
try {Pbn2i) has been observed for most compositions, 
driven by chemical bonding conditions of Ba^+ and R^^ 
cations. The temperature of the P31c -^ Pbn2i phase 
transition was found to vary from 160 K for R=Lu up 
to 355 K for R=l{o.^^^~— Short range magnetic cor- 
relations, evidenced by diffuse neutron scattering, are 
drastically enhanced in the orthorhombic phase, ^^ but 
there is still a controversy about the nature of the low- 
temperature magnetic state. Soda et al.^^ reported only 



short-range correlations in YBaCo407+5 single crystal 
and discuss models where the Kagome and triangular 
layers are fully decoupled, forming spin configurations 
close to 120°. A similar picture of finite magnetic corre- 
lations and decoupled Kagome and triangular spins has 
been adopted to interpret the neutron diffraction mea- 
surements on LuBaCo407+(5 single crystal. ^^ In contrary, 
long range magnetic ordering with an intricate magnetic 
structure has been found for a polycrystalline sample of 
YBaCo407 below 110 K.^ The discrepancy with the sin- 
gle crystal data is possibly due to the oxygen stoichiome- 
try in the crystals studied by Soda et al.^^^^^, an assump- 
tion based on the observation that hyper-stoichiometric 
Y-based and related compositions^ which do not undergo 
the structural transition to the orthorhombic phase, re- 
main short-range ordered. On the other hand, the ^^Fe 
Mossbauer spectra of both YBaC03.94Feo.04O7. 02 and 
YBaC03.94Feo.04O7. 80 polycrystalline samples collected 
at 4K shows well defined magnetic sextets consisting of 
two and three components respectively,^^ a result fully 
compatible with the presence of long-range magnetic or- 
der in both these compositions. 

With the aim to resolve the controversy about the 
nature of the magnetic state in YBaCo407, we under- 
took additional neutron diffraction experiments on sin- 
gle crystal with an oxygen content precisely controlled. 
Our results are consistent with the scenario of long-range 
magnetic order, evidenced by the presence of sharp mag- 
netic Bragg peaks. A model for the magnetic structure 
is proposed from refinements based on a large number 
of independent reflections. In addition, high resolution 
powder X-ray diffraction data clearly indicates that the 
symmetry of the magnetically ordered state is lowered to 
monoclinic P2i. The latter observation testifies of the 
strong spin-lattice coupling in the system and its inher- 
ent link to the magnetic frustration not fully lifted by the 
orthorhombic Pbn2i symmetry. 



II. EXPERIMENTAL PART 

The single crystal of YBaCo407.o was grown by a float- 
ing zone technique in an optical image furnace. A den- 
sified rod of the nominal composition was melted in a 
20% 02/Ar atmosphere, and the crystal was grown at 
1 mm/h. Because of the high affinity for YBaCo407 to 
pick up oxygen on cooling, an afterheater was used to 
keep the growing crystal above 900° C. After growth was 
complete, the atmosphere was changed to 99.995% Ar 
for 24 h and the crystal cooled down to room temper- 
ature over six hours. Oxygen content was verified by 
thermogravimetric analysis of a crushed sample of the 
crystal. The ceramic samples used for the single crystal 
growth and X-ray powder diffraction experiments were 
synthesized from high-purity Y2O3, BaCOa, and C03O4 
reagents by repeated firing in air at 1150°C. The materi- 
als for the X-ray measurements were then treated in ni- 
trogen at 1150°C resulting in the oxygen content 06.95(5) 
as determined by thermogravimetry. 

The neutron diffraction studies were performed on the 
PRISMA time-of-flight spectrometer at the ISIS facility, 
Rutherford Appleton Laboratory (U.K.) and on the four- 
circle diffractometer DIO at the ILL (Grenoble, France). 
For the experiment at ISIS, the single crystal of cylindri- 
cal shape (5 mm diameter, 3 mm height) was mounted 
on an aluminium pin with the [001] vertical axis. For this 
geometry, scattering in the orthogonal (a* 6*) reciprocal 
plane was recorded by rotating the crystal around the 
[001] axis. Data have been corrected for absorption and 
normalized to the incoherent scattering of a vanadium 
standard. The experiment at ILL was done on the same 
crystal prealigned by using the OrientExpress facility. All 
data were collected with an incident neutron wavelength 
of A=2.36 A by using an 80 mm^ two-dimensional mi- 
crostrip detector. Peak integration was performed using 
the program RACER (ILL) in two steps. First, a library 
was built by fitting ellipsoidal shapes to a set of strong 
reflections (/ > Sa), these shapes were used in a second 
pass to integrate all reflections. For each data set the 
list of integrated intensities obtained were corrected for 
Lorentz factors and normalized to the monitor count. To 
refine the magnetic structure (propagation vector k=0), 
432 independent reflections were collected at 5 K and 80 
K. In addition, intensity of some selected peaks was mea- 
sured as a function of temperature (5 K<T<125 K). The 
nuclear structure in the paramagnetic phase was refined 
based on 266 independent reflections measured at 150 K. 
Refinement of both the nuclear and magnetic structures 
was carried out using the FullProf program. ^^ 

A good agreement between the observed and calcu- 
lated structure factors (RF=5.95%) in the paramagnetic 
phase was obtained in the orthorhombic symmetry Pbn2i 
{Pna2i in standard setting) with three types of 120° fer- 
roelastic domains shown in Figured) These domains are 
expected due to the symmetry relation between the or- 
thorhombic Pbn2i and hexagonal P6smc space groups. 
The latter is the common supergroup of Pbn2i and P31c 




FIG. 1: (Color online) A part of the (hkO) reciprocal plane 
demonstrating superposition of the three types of 120° fer- 
roelastic domains, expected in the case of the orthorhombic 
symmetry. 



and is the parent symmetry which has to be used for 
classification of domains and displacive modes in the 
system^ Matrix representation of the rotational parts of 
the three domain generators: {E/000}, {C3+/OOO} and 
{C3-/OOO} were used to introduce the twining law at the 
data processing stage. The low temperature magnetic 
phases at 5 K and 80 K were refined in monoclinic P2i 
space group with six types of domains as discussed in the 
next section. 

High-resolution X-ray powder diffraction data were 
collected on beamline X16C at the National Synchrotron 
Light Source (USA). X-rays of wavelength 0.6911 A 
were selected by a Si(lll) double-crystal monochroma- 
tor. Samples were mounted on a flat brass plate, and 
data were collected in the reflection geometry within the 
temperature interval 10 K<T<200 K. 
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FIG. 2: Integrated intensities of some selected peacks as a 
function of temperature (error bars correspond to 3a interval). 
The indexation relates to the orthorhombic Pbn2i unit cell 
with a -6.30 A, b -10.95 A and c -10.19 A. 
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FIG. 3: (Color online) Example of cj-scan profile for the 
(102) reflection measured at 5 K on the DIO diffractometer 
and fitted with a Gaussian (left) and Lorentzian (right) func- 
tions. The reflection, forbidden in the Pbn2i space group, 
has mainly a magnetic contribution (some small nuclear con- 
tribution can exist since it is allowed in P112i). 



FIG. 4: (Color online) Monoclinic angle as a function of tem- 
perature (left panel). The inset shows the orthorhombic unit- 
cell and the angle 7 which deviates from 90° in the monoclinic 
phase. A part of the X-ray powder diffraction patterns col- 
lected at different temperatures and demonstrating the split- 
ting of the (221) reflection below Ti --100 K (right panel). 



III. RESULTS AND DISCUSSION 

In agreement with the measurements on polycristalline 
materials^ our single crystal neutron diffraction data 
clearly indicate two successive phase transitions in 
YBaCo407, which take place at Ti -100 K and T2 -60 
K respectively (Figl2j); the former has been associated 
with long-range magnetic ordering (propagation vector 
k=0) and the latter with a spin reorientation process.^ In 
the previous single crystal experiments/^^ ^^ the magnetic 
reflections were found to have a Lorentzian shape and a 
width indicating short range spin correlations. The au- 
thors pointed out that the fundamental discrepancy with 
the powder diffraction data^ could be related to a less 
accurate estimation of the magnetic component in the 
powder experiments. In Figure [3l a;-scan for the (102) 
reflection measured at 5 K is presented as an example of 
our single crystal data. This reflection is forbidden by 
the n-glide plane of the Pbn2i space group and there- 
fore has mainly a magnetic contribution. In addition, 
this particular reflection is not affected by the 120° do- 
main structure which can introduce some broadening or 
asymmetry. The reflection has a perfect Gaussian shape 
and a width essentially resolution limited, indicating the 
long-range nature of magnetic ordering. 



A. Symmetry of the magnetically ordered state 

The high-resolution powder X-ray diffraction reveals 
that magnetic ordering is accompanied by structural dis- 
tortions that reduce the crystal symmetry down to mon- 
oclinic (FigH]), producing a clear splitting of some of 
the nuclear reflections below Ti —lOOK. The system- 
atic splitting of the {hkO) reflections and the absence of 



splitting for the {hOl),{Okl) reflections indicate a unit-cell 
with c as unique axis and 7 as monoclinic angle (FiglH 
inset). For simplicity, this non-conventional setting will 
be used to describe the unit-cell rather than transforming 
the coordinates to the conventional cell. The tempera- 
ture dependence of 7 is plotted in Fig. |4](left). Above Ti, 
the small deviation of 7 from 90° is caused by anisotropic 
peak broadening which often precedes magnetostructural 
phase transitions.^^ In fact, above Ti, refinements with 
monoclinic symmetry do not improve what is obtained 
with orthorhombic symmetry, indicating that the struc- 
tural transition coincides with the magnetic one. The 
change of unit-cell metric at the magnetic transition is 
associated with the Cxy component of the strain tensor, 
which spontaneously appears below Ti. This observation 
provides crucial information about the magnetic symme- 
try because a magnetoelastic coupling invariant, linear in 
respect of Cxy^ must be present in the Landau expansion 
of the free energy. Since the bilinear coupling between 
Cxy and the magnetic order parameter, 7^, is forbidden 
by time-reversal, one must consider coupling terms that 
involve even powers of the magnetic order parameter(s). 
Inspection of the irreducible representations (irreps) of 
the Pbn2i space group shows that there are four one- 
dimensional irreps listed in Table [H The strain compo- 
nent Cxy transforms according to the r2 representation. 
It is obvious that only the direct products Fi (g) F2 or 
Fs (g) r4 3 F2 and therefore only coupling invariants of 
the form exyr]ir]2 or exyrjsr]^ are possible. It is alterna- 
tively easy to determine by direct inspection of the irreps 
table that the only way to loose the mirror symmetries 
whilst preserving the two-fold axis along c, necessary to 
stabilize this monoclinic symmetry, is to mix either Fi 
and F2 or F3 and F4. In the former case, the magnetic 
symmetry is P112i and in the latter, PI 12'^. The sys- 
tematic absence of (0 / = 2n + 1) magnetic reflections 



TABLE I: Irreducible representations of Pbn2i space group 
associated with wave vector k=0. 
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(Fig 12]), is uniquely consistent with a time-reversed 2i 
axis, and thus we can conclude that below Ti, the mag- 
netic ordering has the symmetry of the reducible Fs F4 
order parameter. 



B. Model for the magnetic structure 

Both F3 and F4 irreps enters three times in the decom- 
position of the magnetic representation for the 4a sites: 

F(4a) = 3Fi © 3F2 © 3F3 © 3F4 

and therefore one needs to consider three sets of orthogo- 
nal basis vectors for each representation, that are given in 
Table [III The admixture of the two irreps (monoclinic 
symmetry) splits the 4a site into two independent orbits, 
implying that there are no constraints for magnetic mo- 
ments of cobalt sites related by the b and n glide planes 
of the Pbn2i space group, and those can be treated as 
independent variables in the refinement. In contrary, the 
sites related by the 2[ screw axis have parallel in-plane 
{^ixi^iy) configuration and an anti-parallel out-of-plane 
{(fiz) component. 

Another symmetry aspect to consider is the orienta- 
tion domain structure of the magnetically ordered state. 
Since the phase transition involves a change of Laue class, 
mm2 ^2, each domain of the orthorhombic paramag- 
netic phase (FigdJ) gives rise to two ferroelastic domains 

with lattice vectors {a ^b ^ c) and ( a ,- 6 ,- c ), respec- 
tively. Thus, the refinement of the low temperature phase 



TABLE II: Atomic components of basis vectors, (fu 
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Fs and F4 irreducible representations of Pbn2i space group 
entering three times (i=3) in the reducible magnetic repre- 
sentation for an atom occupying the 4a Wyckoff position of 
the Pbn2i space group. 
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FIG. 5: Experimental structure factors vs. calculated ones for 
the refinement of the data collected at different temperatures 
(DIO data). 



must involve six domains and eight nonequivalent Co po- 
sitions. The following strategy for the magnetic struc- 
ture refinement has been applied: initially, only planar 
spin components were allowed to vary (2;-component was 
fixed to zero), in the light of results obtained previously 
on a polycristalline sample^. Given the large number 
of variables, we have imposed the additional constraint, 
not derived from symmetry, that all sites in the Kagome 
sublattice in one hand and in the triangular sublattice 
on the other, have the same moment magnitudes. These 
simplifications, in combination with the symmetry con- 
straints discussed above, leads to 10 independent vari- 
ables. The refinement converged quickly to a spin con- 
figuration equivalent to that proposed from the powder 
diffraction data^^ In the second step, the z-components 
were allowed to vary, with the symmetry constraint dis- 
cussed previously. The refinement included 16 magnetic 
parameters (98 in total) and converged with a good final 
agreement between calculated and measured integrated 
intensities (FigE]) and physically reasonable domain pop- 
ulations. The resulting magnetic structure is shown in 
Figure [6] and the refined parameters are given in Table 
Hm We note that an attempt to refine the magnitude 
of the eight spins independently was not satisfactory, 
most parameters becoming strongly correlated, suggest- 
ing that the best model must be obtained with the afore- 
mentioned restriction on the magnitude of the moments. 
The magnetic model is very similar to that obtained 
from the powder diffraction experiment. The antiferro- 
magnetic ordering in the triangular sublattice is practi- 
cally identical in both models and no z-component for 
these spins was found from the single crystal data. How- 
ever, the spins in the Kagome layers show deviations from 
the planar structure not initially detected in the powder 
measurements, with the largest angle '^40° (FigEl Table 

imi). 

From previous structural work^^ it has been estab- 
lished that the tetrahedra in the Kagome layers make 
a complex tilt pattern in the orthorhombic phase. This 
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FIG. 6: (Color online) Schematic representation of the magnetic structure of YBaCo407 at different temperatures. The spins 
in triangular and Kagome sublattices are presented as large (blue) and small (red) arrows, respectively. 



tetrahedral tilts are not rigid modes of the structure and 
occur with relatively large polyhedral distortions which 
can essentially modify the local anisotropy of the corre- 
sponding Co sites. It appears that the spin configuration 
with a sizable z component for certain sites is due to 
this local anisotropy. The average values of the moments 
in the triangular and Kagome lattices are 2.62(9) and 
1.88(6) \iB respectively, considerably smaller than the 
expected spin contribution of 3.25 /i^ for a Co ion with 
an average +2.25 oxidation state. This suggests that a 
fraction of the moment remains disordered at low tem- 
perature. Additional neutron diffraction data to probe 
a full 2D section of reciprocal space (Figl?]) show that 
indeed diffuse scattering persists down to 2 K. 

The data collected at 80 K have been refined using the 
same strategy. The magnetic structure is shown in Figure 
[6] and the refined parameters are listed in Table Hill The 
model suggests a spin reorientation process which has 
been discussed previously^ based on competing in-plane 
and out-of-plane exchange interactions. The transforma- 
tion from the 5 K structure mainly affects the Kagome 
sublattice and does not change the magnetic PI 12'^ sym- 
metry. Thus, the results for the oxygen stoichiometric 
single crystal are consistent with the powder diffraction 
data obtained for the same composition.^ 

The magnetic structures refined here are obtained with 



less constraints than the powder data and based on 
many independent reflections. They allow to confirm 
that the model proposed in the previous single crystal 
study, ^^ where 120° spin configurations were assumed for 
both Kagome and triangular sublattices is not valid for 
stoechimetric samples. It is likely that the oxygen con- 
tent which can vary largely in these materials is the cause 
of such large discrepancy between experimental results. 
Also, it is important to note that the model presented 
produces a small ferromagnetic moment 0.08(6) \ib per 
Co ion in the (a6) plane. This value is not statisti- 
cally significant but the monoclinic symmetry allows for 
an in-plane ferromagnetic component. It is also diffi- 
cult to determine whether a small ferromagnetic moment 
is present by magnetic measurements which are com- 
plicated by the presence of domains and the disordered 
(glassy) component, the latter producing non-linear mag- 
netization that masks the saturation regime. However 
the results obtained here, allows us to postulate that the 
ferromagnetic moment found in Cai_aji?ajBa(Co/Fe)407 
must arise from a similar non collinear spin arrangement 
rather than from the simple ferrimagnetic ordering with 
antiparallel Kagome and triangular sublattices, proposed 
in Refs^^^^. This assumption is in agreement with the 
recent neutron diffraction study of CaBaCo407 where a 
complex magnetic structure with non-collinear spin or- 



TABLE III: Magnetic structure parameters refined in spherical coordinates (M-value of moment, $-angle witfi tfie x axis and 
B-angle witfi tfie z axis) at different temperatures. Coordinates for Co ions are given as they were refined in the Pbn2i space 
group at 150 K. Domain fractions at T=5 K and 80 K are: 0.34(1), 0.18(1), 0.11(1), 0.10(2), 0.19(2), 0.08(3) and 0.36(1), 
0.18(1), 0.11(1), 0.11(2), 0.16(2), 0.08(3), respectively. 
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"variation of these parameters resulted in their small oscillations, 
within 2°-3° interval near B=90°, indicating that there is no de- 
tectable 2:-component in the triangular lattice and therefore in the 
final refinement these parameters were fixed. 



dering in both triangular and Kagome sublattices has 
been found. ^^ 

Finally, it should be pointed out that the mag- 
netic point group 2' allows a linear magnetoelectric 
coupling and is consistent with a presence of sponta- 
neous electrical polarization along the z-axis. Thus, the 
i?Ba(Co/Fe)407 compositions can demonstrate multifer- 
roic/magnetoelectric properties since the necessary sym- 
metry conditions are fulfilled. 
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FIG. 7: (Color online) Neutron scattering cross sections mea- 
sured in the [ab) plane at different temperatures (PRISMA 
data) . 



Thermodynamics of the phase transitions involving 
coupled order parameters has been discussed in a num- 
ber of works. ^^^^ The relevant case implies two one- 
dimensional order parameters (noted 7^3 and 7^4 to follows 
the labeHng of the irreps) with biquadratic coupling. The 
Landau free energy can be written: 
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(1) 

where ol\^2^ /^i,2 and ^ are coefficients of the expansion. 
The integrity basis of F(r]z^r]^ is formed by the two 
square invariants 7^3, 7^4 and thus, the minimal degree 
to which the expansion (1) has to be truncated to fully 
describe the set of symmetry distinct phases is the fourth 
degree. ^^'^^ However, in this case the transition from the 
high symmetry phase to the phase with 7/3 7^0 and 7/4 7^0 
can only take place in a single point of the phase diagram. 
A more realistic phase diagram, where these phases are 
separated by a first-order transition line, can be obtained 
by including a six-degree term in F (7^3, 7^4) for at least one 
of the order parameters. A comprehensive review of this 
case as well as the case of the symmetric six-degree ex- 
pansion can be found in Refs.^^'^^ However, the isostruc- 
tural phase transition observed in YBaCo407 does not 
appear either asymmetric or symmetric six-degree ex- 
pansion ([1]). Nevertheless, it is straightforward to show 
that eight-degree asymmetric expansion can stabilize two 
distinct thermodynamic phases with 773 7^0 and 774 7^0. 
Minimization of ([1]) in respect of 7/4 gives the equilibrium 



value for this order parameter: 
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Substituting this expression in the F{r]s,r]4) truncated 
at eight degree in respect of r]s gives, with some renor- 
mahzed polynomial coefficients, a well studied case of a 
single one-dimensional order parameter. For this case it 
has been shown^^^^^ that inclusion of invariants with the 
highest degree 2m results in m/2 stable low-temperature 
phases having identical symmetry. Thus, taking into ac- 
count the expression (2), the asymmetric eight-degree 
expansion ^ should split the phase space where both 
components of the reducible Fa F4 order parameter 
are not zero into two distinct regions separated by the 
discontinuous isostructural phase transition. The eight- 
degree expansion necessary to account the isostructural 
transition observed experimentally at T2 indicates that 
the energy difference between the two spin configurations 
is very small. Apparently, this situation is not unique 
for geometrically frustrated systems, where ground state 
degeneracy is lifted by a symmetry lowering. The struc- 
tural distortions can favor several configurations of ini- 
tially infinitely degenerated manifold making the system 
to be easily switchable between the different phases. This 
kind of systems should be sensitive for external pertur- 
bations and may demonstrate rich phase diagrams in the 
coordinates of thermodynamic variables such as pressure, 
temperature and concentration. 



IV. CONCLUSIONS 

An oxygen stoichiometric single crystal of YBaCo407 
exhibits long-range magnetic ordering with propagation 



vector k=0 below Ti ~100 K. The magnetic ordering 
breaks the crystal symmetry Pbn2i to monoclinic 
symmetry P112i which has been observed by high- 
resolution X-ray diffraction. At T2 ~60 K, another 
magnetic phase transition takes place, corresponding 
to a spin re-orientation process and with no additional 
symmetry breaking. In both magnetic phases (above 
and below T2), spins in the triangular sublattice are 
antiferromagnetic aligned in the {ab) plane with no 
detectable z-component. On the contrary, spins in the 
Kagome layers are ordered in a complex non-collinear 
fashion and deviate from the {ab) plane. The transition 
at T2 affects mainly the Kagome sublattice. The average 
values of spins in the triangular and Kagome layers 
are 2.62(9) jj^b and 1.88(6) jub^ respectively at T=5 
K. The reduced values with respect to the expected 
spin contribution for Co in the +2.25 oxidation state is 
related, at least partially, to disorder and short-range 
correlations that persist down to 2 K, as evidenced by 
neutron diffuse scattering. 
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